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We report an extended family of spin textures in coexisting modes of zero-dimensional polari¬ 
ton condensates spatially confined in tunable open microcavity structures. The coupling between 
photon spin and angular momentum, which is enhanced in the open cavity structures, leads to 
new eigenstates of the polariton condensates carrying quantised spin vortices. Depending on the 
strength and anisotropy of the cavity confinement potential and the strength of the spin-orbit cou¬ 
pling, which can be tuned via the excitonic/photonic fractions, the condensate emissions exhibit 
either spin-vortex-like patterns or linear polarization, in good agreement with theoretical modelling. 


Quantized vortices are topological defects occurring 
in many physical systems in optics, condensed matter, 
cosmology and fundamental particles, characterized by 
a phase winding of an integer multiple of 2i r around a 
vortex core. In semiconductor microcavities, quantized 
vortices m and vortex-antivortex pairs Hd may form 
spontaneously in exciton-polariton superfluids and non¬ 
equilibrium polariton Bose-Einstein condensates (BECs). 
Much effort has been devoted to the development of 
methods to create orbital angular momentum(vortices) 
in polariton condensates, providing ways to study the 
fundamental physics of metastable currents or for poten¬ 
tial use as quantum sensors [8] or information encoding 
devices j9]. Optical imprinting m of vortices as well 
as robust spontaneous vortices using chiral polaritonic 
lenses im have been demonstrated. Interestingly, the co¬ 
herent coupling of the photon pseudo-spin (polarization) 
with vortex orbital angular momentum has been shown 
to lead to new types of topological entities, named spin 
vortices, characterised by quantised spin current instead 
of phase winding. Uncontrolled spontaneous spin vortices 
were reported in atomic spinor BECs [12! and in polariton 
condensates subject to structural disorder [13], although 
the exact origin of the polariton spin currents remains 
unclear. We also note that the degrees of freedom asso¬ 
ciated with both the orbital angular momentum and the 
polarization of a photon may find useful applications in 
quantum information processing mm- 

More recently, considerable attention has been focused 
on the investigation of polariton spin-orbit (SO) cou¬ 
pling, i.e. the interaction between the polariton orbital 
motion and its spin due to the effective magnetic field 
induced by the transverse-electric transverse-magnetic 
(TE-TM) splitting characteristic of semiconductor mi¬ 
crocavities m- In condensed matter SO coupling has 
led to significant physical phenomena such as the spin- 
Hall effect [18] and topologically protected conducting 
states [19], whereas in optical microcavities, SO coupling 


of exciton-polaritons enables observations of interesting 
optical counterparts, including the optical spin-Hall ef¬ 
fect [20] . magnetic-monopole-like half-solitons m and 
possibly topological insulators [22j [23]. 

In this paper we demonstrate polariton condensation 
in our recently developed tunable open microcavity sys¬ 
tem [24], where a top concave mirror creates a zero¬ 
dimensional confinement potential for polaritons. Mul¬ 
tiple coexisting condensates are observed under non¬ 
resonant pumping exhibiting an extended family of spin 
vortices and textures. These effects are associated with 
the strong SO coupling in the open cavity system con¬ 
sisting of semiconductor bottom and dielectric top Bragg 
mirrors separated by an air gap. We observe condensate 
emissions showing both spin-vortex-like patterns as well 
as linearly-polarised states. The resultant condensate po¬ 
larization patterns depend on the interplay between the 
strength and the anisotropy of the confinement potential 
and the strength of the SO coupling, which can be mod¬ 
ified with change of exciton/photon fraction. We note 
that polariton condensates exhibiting less rich spin vor¬ 
tex phenomena were observed in a geometry of photonic 
micropillars coupled in a hexagonal pattern [25] . 

The open microcavity system consists of planar bottom 
distributed Bragg reflectors (DBR) and a concave top 
DBR (see Supplementary Information [26] ) controlled in¬ 
dependently by nanopositioners (top-left inset of Fig. [I]), 
which allows free tuning of the spectral resonance by 
changing the mirror separation [24]. A total number 
of 12 GaAs quantum wells (QWs) are grown above the 
surface of the bottom DBR at electric field antinodes, 
allowing the strong exciton-cavity coupling regime to be 
reached with a Rabi splitting of ^ 15 meV [26]. Polariton 
condensation is demonstrated with nonlinear increase of 
emission intensity, sharp linewidth reduction and a small 
blueshift (~ 1 meV) far below the bare cavity mode at 
~ 6 meV to higher energy [26 . 

The top concave mirror induces a strong and al- 
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FIG. 1. Illustrative graph of the new eigenstates formed by 
LGoi mode due to spin-orbit coupling. The TE-TM splitting 
in the cavity lifts the degeneracy of the modes leading to 
three energy levels characterised by spin vortices and anti¬ 
vortices. The mathematical form of the azimuthal part of each 
eigenstate is labeled next to its illustrative diagram. The top- 
left inset is a sketch of the open cavity. The basis of circular 
polarization is defined in the top-right inset. 


most harmonic lateral confinement of the polariton 
condensate [24 . For this reason the system eigenmodes 
are studied in the basis of Laguerre-Gauss modes with 
the SO coupling or the asymmetries in the circular shape 
of the top mirror included perturbatively [26]. In or¬ 
der to fully describe the eigenmodes of the system two 
bases of Laguerre-Gauss modes are needed, one for each 
pseudo-spin component: LG , where cr + /cr _ represent 
polaritons associated with left/right circularly polarized 
light, and p and l are quantum numbers quantifying the 
radial and azimuthal phase evolution, respectively. 

We consider the simplest case of perfectly circular mir¬ 
rors. Since only Laguerre-Gauss modes with l % 0 carry 
non zero orbital angular momentum corresponding 
to a phase rotation of 2 irl either clockwise (phase vortex) 
or anti-clockwise (phase anti-vortex), our analysis starts 
from the first excited manifold (FEM) of the harmonic 
potential LGq±! (without SO coupling): 
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where p s i ( 0 ) is the azimuthal part of the polariton wave- 
function with s = ±1 for polaritons associated with g ± 
polarized light, C(r) is its radial part, and 6 and r are 
angular and radial coordinates, a = y/ Ti/rriLpuJHO > with 
rriLP the lower-polariton mass and ljho the strength of 
the confining harmonic potential. Using degenerate per¬ 
turbation theory and including the SO interaction, one 
obtains the following new eigenmodes (see [26]): 

V’i (r, 0) = 2= C(r ) [ipi-i{0) + ip- ii (9 + tt)] 

fo(r,0) = 2=C(r)[ip n (0) 

V> 3 M) = -Lc(r)[<pu(0) +v?_i_i(6>-7r)] 

+ ( 1 ) 

with eigenenergies: E\ = Eq + a 2 /3tt , = Es = Eo, and 
E 4 = Eq — a 2 / 37 r , where Eq is the energy of the LGqi 
mode and (3 is a parameter describing the strength of the 
SO coupling m- The structure of the new eigenmodes, 
illustrated in Fig. [lj can be understood by observing that 
in the presence of SO coupling the conserved quantity of 
the system is the total angular momentum J = l + s. 
Since l and s are both equal to either 1 or —1 the new 
possible eigenmodes have total angular momentum J = 
—2,0 or +2. The SO coupling lifts the degeneracy by 
coherently combining the J = 0 wavefunctions (pi-i and 
tp-n) to form new eigenstates, while leaving the energy 
of the two J = ±2 modes unaffected. For these J = ±2 
any linear combination of tpn and p- 1-1 is a suitable 
eigenmode in the presence of SO coupling. 

To investigate the properties of spin vortices low tem¬ 
perature photoluminescence (PL) measurements were 
carried out (details in [26]). In the first set of measure¬ 
ments, a concave mirror with a radius of curvature (RoC) 
of 20 gm was employed and the mirror separation was ^ 1 
gm. The cavity is detuned so that polaritons in the FEM 
modes have a photonic fraction of ~ 64%. Below the con¬ 
densation threshold, the spectrum associated with the 
FEM displays two broad features, as shown in Fig. [ 2 ] (a). 
With increase of pump power, condensation occurs and 
the linewidths drop sharply due to an increase of tempo¬ 
ral coherence. Three well-resolved modes labelled by i, ii 
and iii are now revealed in Fig. [ 2 ] (b). Energy resolved 
images, shown in the left panels of Fig. [ 2 ] (c),(d) and (e), 
show a ring-like field distribution for all the three modes. 
The imperfection of the ring shape of mode iii is due to 
slight asymmetry of the confinement potential as will be 
discussed later. A linear polariser and a quarter wave 
plate are inserted into the optical path to collect polar¬ 
ization and energy resolved images for each mode in the 
horizontal-vertical ( 0°/90°) basis, diagonal (±45°) basis 
and circular (cr + /cr _ ) basis, and the associated Stokes 
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parameters, Si, S 2 and S 3 , are calculated for each pixel 
of the image [5B]. The linear polarization angle 0, defined 
as 20 = arctan(S 2 / Si), is mapped out for each mode in 
the middle panels of Fig. [ 2 ] (c),(d) and (e). As the circu¬ 
lar polarization degree (S 3 ) is low for all the three modes 
m, the linear polarization vectors characterize well the 
spin textures. 

All three modes display quantised pseudospin currents 
characterized by a 2 i r rotation of 0 around the mode 
cores, with a high linear polarization degree y/ Sf + S| ^ 
0.95 being exhibited. For both modes i and iii, 0 changes 
nearly linearly with the real space azimuthal angle 6, cor¬ 
responding to the rotation of the vector of linear polariza¬ 
tion clockwise around the mode centre, which indicates 
a co-rotating relation between 6 and 0 , as indicated by 
the right panels of Fig. [ 2 ] (c) and (e). At 6 = 0°, we 
observe 0 =~ 90° (horizontal polarization) for mode i 
and 0 0 ° (vertical polarization) for mode iii ( 0 ° is 

defined as vertical, see the middle panel of Fig. 1(c)), 
showing they are azimuthal and radial spin vortices cor¬ 
responding to the extremal modes 0 i and 04 in Fig. |TJ 
respectively. By contrast, mode ii is a spin anti-vortex 
displaying the opposite pseudospin vector rotation with 
respect to i and iii, with 0 and 6 counter-rotating (right 
panel Fig. [ 2 ] (d)). As discussed in m, its hyperbolic-like 
polarization pattern results from the coherent combina¬ 
tion, with any initial phase difference, of J = ±2 half¬ 
vortices with different polarization (modes 02 and 0 3 in 
Fig-l correspond to the case with a phase difference of 
0 or 7 r). The energy splitting of ^ 0.7 meV observed 
between modes i and iii indicates strong SO interaction, 
consistent with transfer matrix simulations performed for 
the case of a planar open cavity revealing values of TE- 
TM splitting at high momenta of ~ 0.5 meV. Such a large 
value mainly arises from the phase shifts due to reflec¬ 
tions at the air gap interfaces in the open cavity system. 
Possible reasons for the unequal energy spacing between 
modes i, ii and iii are discussed in [26] . 

We also observe spin textures for polaritons condensed 
into higher order LG-associated modes, when these are 
tuned into resonance with the exciton. Similar to the 
LG{)±\ case, SO coupling also mixes modes in the sec¬ 
ond excited manifold like, for example, LG 10 and LG 02 
modes. As illustrated in Fig. [ 3 ] (a), the modes formed 
are quasi-spin vortices labelled as type A and B. The 
polarization vectors exhibit radial (A) or azimuthal (B) 
spin vortex character in the inner core and azimuthal 
(A) or radial (B) spin vortex character in the outer ring, 
connected by transient elliptically polarised states. Such 
quasi-spin vortices of polariton condensates were experi¬ 
mentally observed as shown in Fig. [ 3 ] (b) and (c), with a 
change of linear polarization angle of 7 r between the inner 
core and outer ring. Here above condensation threshold 
four spectrally resolved condensates are observed and for 
simplicity we show polarization patterns only for two of 
them, which fully demonstrate the principle illustrated in 
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FIG. 2. PL spectra of the FEM below threshold (a) and 
above threshold (b) where linewidth narrowing allows the res¬ 
olution of three different modes. The RoC of the top concave 
mirror is 20 pm and the photonic fraction is 64%. (c)-(e) 

Spatial intensity and polarization properties of mode i (c), ii 
(d) and iii (e). The left panels show the real space PL in¬ 
tensity; the middle panels show the real space distribution of 
the linear polarization angle; the right panels shows the lin¬ 
ear polarization angle as a function of the real space winding 
angle circulating clockwise around the white dashed circle in 
the middle panel. 0° and 90° are defined in the middle panel 
of (c) which also defines the spatial coordinates used for all 
figures in this paper. 


Fig. [ 3 ] (a). The imperfection of the mode spatial profile 
and the linear-like polarization vector of the inner core 
in Fig. [ 3 ] (c) compared to (a) are most likely due to the 
slightly elliptical shape of the top concave mirror. 

If the concave top mirror has a sufficiently strong el- 
lipticity which perturbs the harmonic confinement po¬ 
tential along the two orthogonal directions with strength 
a (see [2d for details) it may induce Mathieu-Gauss 
(MG) modes [271128] which are characterised by linearly- 
polarised orthogonal double-lobe profiles (See Fig.Qd)). 
The eigenmodes of the cavity arise from the competition 
between the asymmetry of the mirror and the strength 
of the SO coupling: either spin vortices or linear po- 
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FIG. 3. (a) Schematics of the generation of spin textures 

with higher order LG modes, (b) and (c) Experimental obser¬ 
vation of the non-trivial spin modes in (a), with left, middle 
and right panels showing intensity profiles, linear polariza¬ 
tion angle maps and polarization angle winding. The blue 
and red traces in the right hand panels show the outer and 
inner rotations defined by the dashed white circles in the mid¬ 
dle panel.The horizontal dark lines on Figs 3(b) and (c) arise 
from dark pixels on the ccd detector. 


larised states will be formed depending on which term 
dominates. In order to demonstrate condensation in MG 
modes the strength of SO coupling can be reduced by tun¬ 
ing the energy of the condensed modes closer to the exci- 
ton, and mirrors with smaller RoC chosen where the con¬ 
finement potential is stronger and the spatial anisotropy 
is more pronounced. 

From these considerations a concave mirror with a RoC 
of 7 |itm is chosen, and a photon fraction of 41% employed. 
Fig. [4] (a) shows the spectrum of the polariton condensate 
associated with the FEM, where the low energy modes ii 
and iii are preferentially selected above threshold leading 
to significantly larger intensity than mode i. Non-trivial 
differences, compared to the spin vortices in Fig. [2j are 
found in the mode spatial profiles and polarization pat¬ 
terns, as shown in Fig. [4] (b) and (c) for mode ii and iii. 
Instead of being spin vortices/anti-vortices, modes ii and 
iii clearly show MG-like orthogonal double-lobe profiles 
(left panels) with vertical (ii) and horizontal (iii) linear 
polarization (right panels). The simulated intensity dis¬ 
tribution and polarization maps of one of the eigenstates 
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FIG. 4. (a) Spectrum of the FEM above threshold with 

a photonic fraction of 41% for a top concave mirror with 
RoC=7 pm. (b)(c) Real-space PL intensity profile (left panel) 
and real-space distribution of linear polarization angle derived 
from the Stokes parameters for modes ii and iii (right panel), 
(d) Simulation demonstrating the impact on the spin textures 
of the interplay between SO coupling and the cavity ellipticity. 
The ellipticity parameter a is set to —0.10 meV / fim 2 for all 
graphs. Upper panels show real-space intensity profiles for SO 
coupling parameters (3 = 0.25, 0.05, and 0.01 meV • fim 2 from 
left to right, while the lower panels show the corresponding 
real-space maps of polarization angle. When the SO coupling 
decreases the mode changes gradually from a spin-vortex to 
a linear polarised MG-like profile. The definition of [3 and a 
is detailed in the supplementary materials [26j . 


confined in an elliptical potential are shown in Fig. [4] 
(f) for decreasing TE-TM splitting factors. Theoreti¬ 
cally, it is seen that smaller SO coupling leads to the 
MG mode being the eigenstates of the system as the el¬ 
lipticity term has greater influence. Importantly, for the 
same RoC=7 pm mirror, we can recover the vortex-like 
spin textures for all three modes similar to those shown 
in Fig. [ 2 ] by doubling the photonic fraction up to 82%, as 
shown in the supplementary material [26] . In addition, as 
expected this enhancement of SO coupling with increase 
of the photon fraction results in increase of the i-ii mode 
splitting from 0.71 meV to 1.02 meV. This demonstrates 
the advantages of the tunability of the open cavities in 
permitting the degree of the light/matter fractions of the 
cavity polaritons to be varied but also in allowing flexible 
manipulation of the condensate polarization textures. 

In summary, we have demonstrated polariton conden¬ 
sate emission exhibiting spin vortices and more elaborate 
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spin textures in a tunable microcavity system with lateral 
confinement. We note that below threshold the tomo¬ 
graphic energy-resolved images of broad polariton modes 
also exhibit spin-vortex structures, but not all patterns 
are observed due to the lack of spectral resolution. In 
order to reveal the complete mode structure it is im¬ 
portant to be in the regime of condensation where the 
temporal coherence is increased. Our observations are 
promising for optical information applications involving 
photon orbital angular momentum and polarization [9] 
and are interesting for further studies on cavity quantum 
electrodynamics (CQED) involving polaritonic spin. 
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FIG. SI. Schematic of the experimental setup. 


SAMPLE PREPARATION AND EXPERIMENTAL SET-UP 


The open microcavity system used to perform the experiments consists of a 31-pair Alo. 2 Gao. 8 As/Alo. 95 Gao. 05 As 
bottom distributed Bragg reflector (DBR) with a near-surface active region and an 11-pair Si 02 /Ti 02 circular-shaped 
concave top DBR separated by a micrometer sized gap. Arrays of concave mirrors are fabricated through focused 
ion beam (FIB) milling of a planar Si 02 substrate before coating with dielectric layers [29;. The radii of curvature 
of the concave mirrors used in this work were 20 qm and 7 qm and can confine the polariton mode down to 1-2 
qm. Nanopositioners allow independent positioning of both DBRs to form a planar-concave cavity where the spectral 
resonance can be tuned by changing the separation between them (Fig. |S1[ left inset) [24] . In the cavities used for this 
experiment, three sets of four 7 nm GaAs quantum wells (QWs) are embedded in the active region at the antinode of 
the optical field to enhance the exciton-photon coupling. 

The experimental setup is shown in Fig. [ST] The cavity system is placed in a vacuum tube with a small amount 
of He exchange gas which is immersed in a liquid helium dewar. Optical access to the cavity is provided by placing 
an optical table on top of the dewar. The sample is non-resonantly excited at 630 nm, close to a stopband minimum, 
with a spot size of ^ 30 qm on the top mirror surface. The beam is reflected into the dewar by a beam splitter (BS) 
to an objective lens above the sample (NA=0.55) and the photoluminescence from the cavity is collected along the 
same optical path. The final image is sent to the end facet of a wound fibre bundle (WF) consisting of a 4x4 mm array 
of single mode fibres, which is imaged onto the spectrometer slits. It should be noted that light is depolarized by the 
WF, avoiding any measurement error possibly induced by the polarization-associated efficiency of the spectrometer 
mirror/gratings. Fig. SI represents the setup used for k-space imaging, where lenses are located both in the dewar 
and on the optical table to form a confocal imaging system which projects the k-space image onto the WF facet. This 
setup can be easily altered for real-space imaging by simply replacing the two lenses on the optical table by one that 
focuses the real space image on to the WF facet, as used for our studies on spin vortices. A linear polarizer, whose 
polarization axis can be varied by motor-controlled rotation, is inserted before the WF to obtain polarization-resolved 
images, while adding a A/4 plate enables acquisition in circular polarization basis. 
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FIG. S2. Spectra of LPB for 0D cavity as a function of piezo voltage on bottom DBR. The simulated bare exciton (X), bare 
cavity (Cav), LPB and UPB dispersions are labelled, showing a Rabi splitting of 15 meV. 


DEMONSTRATION OF STRONG COUPLING 

The strong coupling of the 0-dimensional (0D) cavity is demonstrated by scanning the PL spectrum with varying 
cavity length, as shown in Fig. |S2| The cavity length is decreased by applying a DC voltage to the bottom z- 
nanopositioner, raising the bottom sample closer to the top mirror. The change in cavity length as a function of 
applied piezo voltage is closely linear for voltages less than V = 20 V. The lower polariton branch (LPB) shows 
strong curvature as a function of voltage arising from anti-crossing with the exciton energy, the signature of strong 
exciton-cavity coupling. The upper polariton branch (UPB) is not observed due to the strong absorption induced by 
the excitonic continuum of the 12 GaAs quantum wells. Fitting the LPB dispersion curve of the LGoi mode with 
a coupled oscillator model gives a Rabi splitting of ~ 15meV, in good agreement with theoretical and experimental 
values reported in similar structures mum]. 

The Rabi splitting can also be estimated from the 2-dimensional polariton dispersion. As shown in Fig. 1(a) 
of Reference [24] the top DBR contains planar regions around the arrays of concave features. The formation of a 
planar-planar cavity using this region gives rise to two-dimensional polaritons. Fig. |S3| (a) and (b) shows the angular 
dispersion of the planar open cavity at two-different detunings of —9.4 meV and —1.5 meV. In both cases the LPB 
can be fitted with the expected Rabi splitting of 15 meV, in agreement with previous full microcavities containing 
similar numbers of QWs mm- 
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FIG. S3. Angular resolved spectra of LPB with planar top DBR for detunings of —9.4 meV (a) and —1.5 meV (b). In both 
cases the simulated dispersions show a Rabi splitting of 15 meV. 


(a) (b) (c) 



FIG. S4. Power dependence of the first excited manifold (FEM) of states at S = —4 meV, which splits into a triplet above 
threshold (black, red, blue traces), with spectrum in Fig.2 (b) of the main text. Integrated intensity (a), spectral linewidth (b) 
and energy blueshift (c) are shown for all the three modes. 


POLARITON CONDENSATION 


Polariton condensation into the first excited manifold (FEM) of states is characterized in Fig. S4 Non-linearity of 
output intensity and sharp linewidth reduction is observed. The blueshift at threshold, ~ 0.8 meV, is significantly less 
than the E cavity — Elpb which is ^ 6 meV, showing that the cavity is in the strong coupling regime above threshold. 
The spectra below and above threshold are in Fig. 2 (a) and (b) of the main text. 
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FIG. S5. Tomographic images of modes i, ii and iii in Fig. 2 of the main text, unpolarized (unpol), and with a linear polarizer 
with the polarization axis oriented at 0°, 45°, 90° and 135°. 


EXPERIMENTAL DETERMINATION OF THE STOKES PARAMETERS 


In our experiment polarization-resolved images are taken in order to obtain the spatial intensity distribution in the 
horizontal-vertical (0° and 90°), diagonal (45° and 135°) and circular (a + /a~) bases. Fig. S5 shows the clear rotation 
of the intensity distribution as a function of the linear polarizer angle. Modes i and iii show orthogonally positioned 
lobes co-rotating with the polarizer axis while mode ii anti-rotates with the polarizer axis. This indicates the formation 
of spin vortices for modes i and iii and spin anti-vortices for mode ii. Stokes parameters for horizontal-vertical (Si), 
diagonal (S2) and circular (S3) basis are obtained for each spatial pixel of the mode by 


1(0°) - J(90o) 
1(0°) +7(90°) 
1(45°) - /(—45°) 
7(45°) + J(-45°) 
J(q + ) - J(cr-) 

I(<7+)+/(<7-) 


where 1(0) is the measured intensity in various polarization basis. The polarization angle in real space 0, is obtained 
from 20 = arctan(S2/Si), as the winding angle in the S1-S2 plane of the Stokes presentation is twice of that in real 
space. Fig. |S6| shows the Stokes parameters derived for the spin vortices in Fig. 2 of the main text. 
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FIG. S6. Stokes parameters S'i (a), S 2 (b) and Senear = \JS\ + S'| (c) derived for Modes i, ii and iii in Fig. 2 of the main 
text. 


SPIN VORTICES WITH A TOP MIRROR OF ROC=7 pM 


The mode spatial profiles and polarization patterns with the RoC=7 pm top mirror with a high photonic fraction 
(82% photon-like) are shown in Fig. S7 Azimuthal, radial and hyperbolic spin vortices are revealed, similar to Fig. 2 
but contrasting sharply with Fig. 4 of the main text taken at 41% photon fraction where the ellipticity perturbation 
plays a larger role relative to the SO interaction. 
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FIG. S7. Spatial intensity and polarization properties of the three eigen modes with a RoC=7 pm top mirror, shown in (a),(b) 
and (c) respectively for a photonic fraction of 82%. The left panels show the real space PL intensity; the middle panels show 
the real space distribution of the linear polarization angle; the right panels show the linear polarization angle as a function of 
the real space winding angle circulating clockwise around the white dashed circle in the middle panel. 


THEORETICAL MODEL 


We present here the theoretical approach used to interpret the experimental data. In our approach we use degenerate 
perturbation theory to find the eigenmodes of the polariton system in the case of low polariton densities (i.e. when 
the nonlinearities play a negligible role) in the presence of SO coupling, elliptical shape of the top concave mirror and 

birefringence ten fhe anisotropy of the refractive index of the top mirror. Taking the eigenvectors (J) and (J 

to represent cr + /cr _ circularly polarised polaritons, the 2x2 Hamiltonian describing the lower-polariton branch in the 
linear regime can be written as: 


H = 


_ jrv_ 

2rriLF 


■V 


*<!)' +^"72 -££ + '■' 


(SI) 


where tulp is the lower-polariton effective mass. The terms depending on (3 = ?i 2 (1 / rn^p — l/mJIO/2, where rrt^p^™ 
are the lower-polariton masses in the TE/TM polarizations, describe the TE-TM splitting. As pointed out in [24] the 
top concave mirror induces a strong near-harmonic lateral confinement potential. V = \rriLpu 2 HO [x 2 (l+5)+y 2 (l — £)], 
where ujho is the strength of the harmonic confinement. The terms +a = +\mppw 2 H0 5 account for an elliptical 
asymmetry of the top circular mirror with the long axis either aligned along the x or y directions. As birefringence 
may arise in both the top and bottom mirrors due to strain, the terms Qe ±l6 /2 account for a birefringence that 
induces a shift at k = 0 between the TE-TM branches and tends to align the field polarization along the direction 6 . 

Since the harmonic confinement is much stronger than the SO coupling, the birefringence and the asymmetry, one 
can treat these terms as perturbations. To study a 2-dimensional harmonic oscillator several equivalent eigenvector 
bases can be used. Among them two are particularly useful: the basis of Laguerre-Gauss modes LG ^ (where p 

and l are radial and azimuthal quantum numbers) and the basis of Hermite-Gauss modes HG°± (where r and s are 
quantum numbers along the x and y axes). While the basis of LG modes allows a more intuitive understanding of 
the shape of the spin vortices, the basis of the HG modes allows an easier evaluation of the matrix elements needed 














to determine the perturbed eigenenergies and eigenmodes. For this reason, and since the perturbed eigenmodes and 
eigenenergies do not depend on the basis of the Hilbert space used to evaluate them, we use the basis of the HG 
modes to apply perturbation theory. In the case of the first excited manifold the four relevant HG modes are: 


ipi (x,y) 
•03 (x,y) 


V<rV2 

* 2 +y 2 

ye 2 <j 2 
y <t 4 7t/2 



v >2 (x,y) 
V>4 (x,y) 


xe 2 cr 2 
\J (j 4 7r/2 


ye ^ 
y<7 4 7r/2 



Using these modes as basis, the new perturbed eigenenergies and eigenmodes of the system are obtained by diago- 
nalising the following matrix: 




M = 


fairer 6 


|7ra 2 (-2/? + e-*V 2 U) 

0 

— §7 v/3a 2 


7T(T 2 {-2f3 + e ld (j 2 n) 


|a7rcr 6 

§7 T/3a 2 


0 

— |7t/3ct 2 


— fairer 6 
\iuj 2 (2f3 + e~ ie a 2 Q) 


17T^(J 2 

0 

\tt<j 2 (2(3 + e 2< V 2 U) 


— ^a7rcr 6 


(S2) 


where cr = ^Ti/itllpujho- The eigenmodes and eigenenergies of this matrix reduce to those of Equation (1) in the 
manuscript in the case of zero birefringence and no asymmetry in the harmonic confinement. For this particular case 
the energy spectra, the polariton density and the polarization angle for the four eigenmodes are plotted in figure 
S8 As expected, the higher (black) and the lower (red) modes are azimuthal and radial spin-vortices respectively, 


while the two remaining central modes (green and blue) are spin-antivortices in agreement with the experimental 
observations in fig 2 (c,d,e) in the main text. 
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FIG. S8. Case for no mirror ellipticity and birefringence equal to zero, (a) Simulated energy levels: mode i (black), mode 
iia (blue), mode iib (green), and mode iii (red), the energy is evaluated with respect to the unperturbed mode, (b) Simulated 
polariton density (first row) and angle of linear polarization ardan{S^l Si) (second row). To obtain both the spectra and the 
eigenmodes the following parameters are used: /3 = 0.2 meV ■ /im 2 , a m 0.7 ym. All graphs are 4 ym x 4 ym in size. 

The case of asymmetry and birefringence different from zero is plotted in Fig.|S9[ Clearly, the effect of these terms is 
to lift the degeneracy among the two central modes, thus breaking the symmetry of the spectra and inducing polaritons 
to polarise along a preferred direction. This is consistent with what it is observed in the experiments, although in the 
experiments the shape of the high-energy mode is generally less deformed by the asymmetry and birefringence than 
the low-energy modes. A possible explanation for this is that our theoretical model is based on the approximation 
of quadratic dispersion while in the polariton system the dispersion is strongly dependent on k. Since the modes 
are strongly confined, high k vectors are likely to play an important role. In addition exciton-exciton interactions 
and pump-decay mechanisms, both of which are not included in our model, may also lead to experimental/theory 
differences in the details of the patterns. 
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FIG. S9. Case with mirror ellipticity and birefringence different from to zero, (a) Simulated energy levels: mode i (black), 
mode iia (blue), mode iib (green), and mode iii (red), the energy is evaluated with respect to the unperturbed mode, (b) 
Simulated polariton density (first row) and angle of linear polarization arctan(S 2 / Si) (second row). The following parameters 
are used: — 0.2 meV • /im 2 , a = 0.7 /im, Q = 0.4 meV, 6 = 0.017T, and a = —0.6 meV to obtain both the spectra and the 
eigenmodes. All graphs are 4 fim x 4 /im in size. 
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Polarization Map 


FIG. S10. Simula ted characters of the quasi-spin vortex A of Fig. 3 of the main text, including amplitude, linear polarization 
degree \f S* + circular polarization degree S 3 and linear polarization angle arctan(S 2 /Si) /2. These characters illustrate 
the corresponding polarization pattern presented in the main text. 


SPIN TEXTURES GENERATED FROM HIGHER ORDER LG MODES 

We derive the polarization patterns by superposing cr + LGio and a~ LG 02 • The result for the type A quasi-spin 
vortex in Fig.3 of the main text is shown in Fig. |S10| as an example. The polarization pattern shows a radial spin 
vortex in the inner core and an azimuthal spin vortex in the outer ring, in good qualitative agreement with the 
experimental data shown in Fig. 3 (b) of the main text. 



















